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Non-Perturbative Renormalization and the Fermilab Action * 

Huey-Wen Lin a 
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We discuss the application of the regularization independent (RI) scheme of Rome/Southampton to determine 
the normalization of heavy quark operators non-perturbatively using the Fermilab action. 

1. Introduction 



Heavy meson hadronic matrix elements play 
an essential role in determining many of the pa- 
rameters of the Standard Model, such as the 
CKM matrix elements. Theoretically, lattice 
gauge theory provides a means of computing 
these hadronic matrix elements with better con- 
trol of systematic errors yQ. The heavy quark 
methods that have been successfully used on the 
lattice are: Shcikholcslami-Wohlert (SW) action, 
nonrclativistic QCD, heavy quark effective theory 
(HQET) and Fermilab action. Here we focus on 
the Fermilab action because it may best control 
the systematic errors associated with the large 
quark mass. 

Renormalization of lattice operators is neces- 
sary in order to obtain physical results from nu- 
merical simulations. Non-perturbative renormal- 
ization (NPR) methods are attractive since they 
avoid possibly large errors associated with lattice 
perturbation theory. This is especially important 
when domain wall fermions are used for the light 
quarks. Among various NPR methods, we decide 
to focus on RI/MOM scheme NPR g|, because 
of the promising results [3| reported by the RBC 
collaboration. We will use NPR to represent the 
RI/MOM scheme NPR method in the rest of this 
paper. 

Our goal is to adopt the RI/MOM scheme and 
to extend its application to heavy quarks. In par- 
ticular we will show how to generalize the usual 
NPR "window restriction" : Aq CD << — p 2 « 
1 /a 2 to the kinematic region required by the Fer- 
milab approach to heavy quarks. 



P2- q 




Figure 1 . Example of a Feynman diagram in weak 
interaction effective field theory. 

2. Generalizing RI NPR to heavy quarks 

We need to find a new renormalization re- 
gion, which meets three requirements: i) It per- 
mits evaluation in continuum perturbation the- 
ory. Thus, it should involve non-exception exter- 
nal momenta, giving no infrared singularities as 
light quark and gluon masses are set to zero, ii) 
The Fermilab method remains valid, i.e. the re- 
gion is not too far from the heavy quark mass 
pole, iii) It can be studied using a Euclidean, 
lattice Green's function. 

Let's take a look at a typical Feynman diagram 
with a one loop QCD correction integral, for ex- 
ample, that shown in Figure ^ with the integral 



l q 1 i(#i-j^ + m) i(v/2-i + m) 
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q) 2 - m 2 (p 2 - q) 7 



(1) 



Here q is the internal gluon momentum and pi 
is the net external quark momentum flowing into 



internal quark line within the loop. 



"Talk presented in Lattice 2003. 



The denominator of the i quark line, with 
Minkowski external momenta but Euclidean loop 
momentum is (pi—q) 2 — p 2 + q 2 + m 2 — 2ipifl-qo. 
Since qo may be small, it doesn't guarantee that 



2 



Space-like 



A 2 QCD « « 1/a 1 



m 

Time like 



Figure 2. The complex-p 2 = p$ — p 2 plane. 



we can avoid the pole on the po axis during inte- 
gration. Therefore, additional conditions have to 
be set to avoid singularities and permit the use of 
continuum perturbation theory. Specifically, we 
require that the absolute value of each denomi- 
nator be much larger than Aq CD . 

Let's first check our formalism with a massless 
quark. We find that — p 2 = [i 2 » Aq CD , which 
is consistent with the the original NPR method 
for light quarks. 

How about heavy quarks? We propose: m 2 — 
Pq = /i 2 >> A.q CD . In order to keep Fermilab 
discretization errors under control, we have to 
constrain the above quantity to be much smaller 
than 1/a 2 as well so that the heavy quark lines 
remain nearly on-shcll. Therefore, our "window" 



of rcnormalization is A 2 - 



Figure [2] is a plot in complex p 



qcd «m 2 -p 2 « 1/a 2 . 

J plane. 



Our choice of renormalization points are located 
within the small region circled. The distance from 
this circled region to the mass pole has to be much 
smaller than and much larger A.q CD for per- 
turbation theory to be valid. This is what I mean 
by "slightly off-shell condition" in the rest of the 
paper. 

Note that the region we choose is time-like, 
while the usual NPR's choice is space-like. Why 
do we pick a time-like point, instead of a space- 
like point? The constraint comes from the Fer- 
milab action, which is an on-shell 0(a) improve- 
ment. If we choose a deep, space-like point, the 
well-behaved properties of Fermilab action would 
no longer be valid. 



3. Off-shell Fermilab action 

Recall "heavy quark" means the product of 
m and lattice spacing a may be larger than or 
equal to 1. Under this circumstance, it's natural 
to expect an improved action to break the time- 
space symmetry. Here is the original Fermilab 
action @|: 

S = ^ i>ni>n 
n 

-Kt^2bij„(l - 7o)t r n,oV'„ + 6 + ^n+ai 1 + 7o)Uj,, W 

n 

~K s ^2[^ n (r s - 7i)U Tl>j + $ n+ i(r s + liWn^n] 
n , i 

+ic E K s ^ ^Pn^OiEn-i^n (2) 



Note that those coefficients are the functions of 
ma, and are expected to remain finite for both 
small and large ma. 

What happens when we go slightly off-shell? If 
the external quark masses differ from the physi- 
cal ones by a small amount 5m, with aSm <C 1 
then the resulting errors will be no larger than 
the other discretization errors in the Fermilab ap- 
proach. 

However, when we go off-shell, we have to in- 
clude more terms to make the 0(a) improvement 
complete. (Why? because contact terms and new 
non-gauge-invariant operators can appear.) First, 
we have to add 0(a) off-shell improvement terms 
in the action. Fortunately, those terms which 
need to be added in the action can be compen- 
sated by improving the quark fields 5 . There- 
fore, the action itself remains the same, including 
the well-behaved coefficients. 

The improved quark fields will look similar to 
those in ref [5], but with broken space-time sym- 
metry. Likewise for the composite operators. 
The broken symmetry gives us at least double 
the number coefficients to be determined non- 
perturbatively, compared with the light quark 
cases. We are still exploring how these coefficients 
may be determined in a more efficient way. 
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4. Imposing RI NPR on the lattice 

To respect the "slightly off-shell condition" , we 
need to change the time-component Fourier trans- 
formation to Laplace transformation. Take the 
operator Or(x) — tp(x)Ttp(x) as an example, for 
a general Dirac matrix T. The non-amputated 
Green's function in momentum space is: 



G ( P a) = {S {U} ( P | 0) T j 5 S\ u} (p | 0)75 



\{U} 



(3) 



where Si(p \y) = Jd 3 xJ dx Si{x | y)e +l P- x e PoXo . 

We calculate the amputated Green's function 
and apply the projector onto P Q to define: 

IoM = ^Tr[S(pa)- l Go{pa)S{pa)- l P ], (4) 

where the S(pa) is the gauge-averaged propaga- 
tor. Following the usual NPR procedure, we ap- 
ply the renormalization condition: 

{p\Or R \p) U_ p2= ^ = (p\Or\p)o, (5) 

where (...) represents tree-level value. Thus, im- 
posing ZoZ^' 1 ro| m 2_ p 2 =M 2 = 1 will determine 
the renormalization factor Zo{na,g{a)). 

5. Matching to the MS scheme 

Our lattice NPR calculations are done in the RI 
renormalization scheme. We need continuum RI 
renormalization calculations to match to physical 
values given in the continuum MS scheme. The 
previous matching calculations done at — p 2 = /i 2 , 
while taking m — > must be extended to our new 
scheme. 

Let's take a look at the simplest example: the 
calculation of Z q and Z m . Our new renormaliza- 
tion conditions become: 



i dm.jp)) 



= i 



— Z q Z m Tr[Y, 2 {pj\ m 2- p 2 =l _ l 2 = 1 



(6) 



(7) 



where i5' _1 (p) = Z q \j>T. 1 {p) - Z m mT. 2 {p)}. 

The calculations are done in the dimensional 
rcgularization scheme, where D is 4-2e. Here are 
the matching factors for Z m and Z q in the RI and 



MS schemes: 

Z™ a s N 2 -l ( 8 + 3£ m 



Z MS 



4tt 2N c 



m 2 — /i 2 



2 4 2 

-[( 3 + 0-5—2 - ^7^2—2^]^" 



m 2 — (m 2 — [i 2 ) 2 fi 



,' 2 



-Hn^(3 + 20} + O(ai 



(8) 



yMS 

q 



4tt 2N c sv 2 to 2 - ^ 2 



m 2 , m' 2 



where // is defined a,s g = Z g g[i' e . 

6. Conclusion 

In this paper, we present a first NPR renormal- 
ization proposal for relativistic heavy quarks. We 
propose a new renormalization point: A 



QCD 



« 



m 2 — p 2 « I /a 2 , in a slightly off-shell region. 
We must add off-shell improvement terms and de- 
termine the corresponding coefficients as well as 
evaluating Laplace transformations in the time. 

The author would like to thank each member of 
the RBC collaboration for their assistance. Spe- 
cial thanks go to N. H. Christ for inspiration and 
constant discussions on various topics. 
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